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1 Introduction 

In this paper we present an asymptotic formula for the eigenvalue counting function 
of the Schrodinger operator —A + V for unbounded potentials V on several types of 
unbounded fractal spaces. Such an asymptotic formula is often attributed to Niels Bohr in 
the Euclidean setting. We identify a set of sufficient conditions for Bohr’s formula to hold 
on locally self-similar metric measure spaces which admit a cellular decomposition, and 
then verify these conditions for fractafolds [55,60] and fractal fields [25] based on nested 
fractals. In particular, we are able to partially answer a question of Fan, Khandker, and 
Strichartz [19] regarding the spectral asymptotics of the harmonic oscillator potential on 
the infinite blow-up of a Sierpinski gasket (abbreviated SG). 

All these results have similarities in the classical theory of ID Sturm-Liouville operators 
(see [ 18]). The deep analogy between nested fractals (the typical representative being SG) 
and the real line Mi. = [0, oo) is related to the fact that all of them are finitely ramified. 
(A set is said to be finitely ramified if it can be divided into several disconnected subsets 
upon removing a finite number of points from the set. For it suffices to remove one 
point; for SG, two points.) 

Let us recall several known results from the spectral theory of the ID Schodinger 
operator 

Hip = —ip" + V(x)ip, x>0 (1.1) 

with boundary condition at x = 0 of either Dirichlet type, ip( 0) = 0, or Neumann type, 
ip’( 0) = 0. 
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I. Assume that V(x) —> +oo as |x| —> +oo. Then, by the result of H. Weyl, the 
spectrum of H in L 2 ([0, oo), dx) is discrete and, under some technical conditions, 

1 r°° 

N(\,V) := #{A i(H) < A} ~ - / y/(X -V(x)) + dx. (1.2) 

^ Jo 

This is known as N. Bohr’s formula, see [28,34,36]. 

II. Assume that V(x) is compactly supported, or (weaker assumption) vanishing fast 
enough (see below). Put V{x) = V + (x) — V-(x), where V + = max(0, V) and V- = 
max(0, —V), and 

1V_(R) :=#{A; <0}<1V_(-F_(-))- (1.3) 

The estimate of N_(V) as a result can be reduced to the negative potentials (po¬ 
tential wells). We use the notation N-(V) assuming here that V(x) = —V-(x) < 0. 
The following estimates of N-(V) are popular in applications (see [48]): 

(a) (Bargmann) 

roo 

N-(V) < 1 + / xV{x)dx. (1.4) 

Jo 

(b) (Calogero) If V(x) decreases with |x| as |x| — > oo, then 

roo 

N-(V) < Co / y/V(x) dx. (1.5) 

Jo 

The Calogero estimate has the correct scaling in the following sense. 

(c) Consider the operator 

Hfj'ip = — ip" T uVq(x)' 0, x > 0 (plus boundary condition). (1.6) 

Then as a —> oo, 

roo 

N-(crVo) ~ cicr 1 / 2 / y/V(x) dx. (1.7) 

Jo 


This is the so-called quasiclassical asymptotics. It is an important problem to 
find such an estimate for IV_(V), which has true scaling in M. d , d > 2, i.e., for 
any a, 


N-(aVo) < a d/2 $(V 0 ). 

(Cwickl-Lieb-Rosenblum) 

(1.8) 

For d > 3 this is the CLR estimate 



N-(V)<c d I 

\V(x)\ d/2 dx. 

(1.9) 

m d 




For d = 2 the recent results by Grigor’yan and Nadirashivili [22] and Shar- 
gorodsky [53] give the desirable (though not simple) estimate. The paper [ 13] 
contains the justification of the physical conjecture by Madau and Wu on 1V_( V) 
for 2D operators. The case d = 1 was studied in the relatively recent papers 
by K. Naimark, G. Rozenblum, M. Solomyak et al (see [37, 52] and references 
therein). 
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Figure 1: Part of an infinite blow-up of SG( 2), which is Type (i) of the fractafold considered 
in 56.2. 



Figure 2: The ladder periodic fractafold based on SG( 2), which is Type (ii) of the fractafold 
considered in 56.2. 


In this paper we address the item I. above in detail. Items II. a, II. b, II.c will be the 
subject of future work. 

Our main objective is to consider, instead of the Euclidean space, a fractafold, which 
according to Strichartz [ ] is defined as “a space that is locally modeled on a specified 

fractal, the fractal equivalent of a manifold.” The first instance of a fractafold is the 
infinite Sierpinski gasket (Figure 1). As shown by Barlow and Perkins [8], the heat kernel 
Pt(x,y) on the infinite Sierpinski gasket satisfies a sub-Gaussian estimate with respect to 
the Euclidean metric d(-, •): 


p t {x,y)-cit ds/2 


exp — C 2 


d(x,y)° 


d„, \ V 


where d s = 2 log 3/ log 5 and d w = log 5/ log 2 > 2. Here x means that there are upper and 
lower estimates, but the constants ci, C 2 in them may be different. We would like to note, 
however, that the heat kernel is not immediately relevant for spectral analysis, partially 
because its form is complicated, but mostly because the domain of the Laplacian is not an 
algebra under multiplication [11]. Other typical examples of fractafolds that we consider, 
see [60] and Section 6 for details, are shown in Figures 2, 3, and 4. For background 
concerning spectral analysis on fractafolds, see [2,4,15,19,29-31,43-45,47,49-51,56,61]. 
Existence of gaps in the spectrum is investigated in [9,17,2 , ]. Wave equation on 

fractals is discussed in [3,13,14,41,59]. Physics applications, and spectral zeta functions 
in particular, are given in [1,10,16,18,27,54,62], 


2 Main results 

2.1 Spectral asymptotics of — A + V 

In all the examples to follow, K is a compact set in M, d endowed with a Borel probability 
measure p and a “well-defined boundary” dK which has //-measure zero. We shall assume 
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Figure 3: The hexagonal periodic fractafold based on SG( 2), which is Type (ii) of the fractafold 
considered in 56.2. 



Figure 4: The triangular lattice finitely ramified fractal field based on SG( 2), which is considered 
in §6.3. 

that there exists a well-defined self-adjoint Laplacian operator —A A (resp. —A v ) on 
L 2 (K, h) satisfying the Dirichlet (resp. Neumann) condition on dK. Note that dK might 
not coincide with the boundary of K in the topological sense. We assume (as is well 
known in examples) that both —A A and —A v have compact resolvents, and hence have 
pure point spectra. It then makes sense to introduce the eigenvalue counting function 

N h (K , /j. A) := # {A,(-A b ) < a} , b G {A, V}. (2.1) 

Assumption 2.1. There exists a positive constant d s such that 

0 < lim X~ da ^ 2 N h (K, /j,, A) < lim X~ da/2 N h (K, p, A) < oo, (2.2) 

A—>-oc A ^-°° 


where b G {A, V}. 

A stronger condition than Assumption 2.1 is 

Assumption 2.2 (Weyl asymptotics of the bare Laplacian). There exist a positive con¬ 
stant d s and a right-continuous with left limits (cadlag), T-periodic function G : M —> M+ 
satisfying 

(Gl) 0 < inf G < sup G < oo. 

(G2) G is independent of the boundary condition b G {A, V} 
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such that as A —> oo, 


N h (K,p, A) 


\ds/2 


G 


-logA) +A b (A) 


(2.3) 


where R h ( A) denotes the remainder term of order o(l). 

Remark 2.3. The parameter d s is often identified with the spectral dimension of the bare 
Laplacian —A on L 2 (K,p). If K is a domain in with a nice boundary, and p is the 
Lebesgue measure, then d s = d and G is an explicit constant (2ir)~ d p(B)p(K) , where B 
is the unit ball in M d . However, there are classes of fractals K for which (2.3) holds with 
G being possibly nonconstant. 

In many examples, the leading-order term in A b (A) gives information about the bound¬ 
ary of the domain. For an Euclidean domain in with nice boundary, the leading-order 
term of A b (A) scales with A^ 1 / 2 , and the sign of this term is negative (resp. positive) if 
b = A (resp. if b = V) [L2,21,32,35], 

We now consider an unbounded space AAx, which admits a cellular decomposition into 
copies of K. Formally, let K^ := U a K a , where 

• Each K a is isometric to K via the map (f> a : I\ —> K a . 

• We identify dK a := (j> a (dK) to be the boundary of K a , and A° := K a \8K a the 
interior of K a . 

• (Cells adjoin only on the boundary.) For all a ^ a! , (K a n K a i ) = ( dK a H dK a i ). 

Let Ha := p o (j)~ l be the push-forward measure of /j onto K a . For any a ^ a' , it is direct 
to define the “glued” measure pL Q ,a' on Ka U K a i in the natural way: 


MB G B(K a U K a r) : /J- a ,a'(B) = n a (B n K a ) + ii a '(B n K a >). (2.4) 


By extension we define the measure on K 00 . 

Proposition 2.4 (Decoupling of L 2 ). For all a / ol we have K° x fl I\°, = 0 and L 2 (J\° U 
K°,,Pa,a') = L 2 { K aiha) © L 2 (K°,, p a >). 


Proposition 2.4 allows one to decouple the Laplacian on the glued measure space 
into the direct sum of the Laplacians on the individual components (see [ 18, Proposition 
XIII.15.3]): 


A b -A b 


A b 


(2.5) 


from which it follows that 


N h (K a U A) = N h (K a , p ai X) + N h (K a i, p a >, A). (2.6) 

By extension we have that 

N h (K 00 ,p 00 ,X) = J2 Nh ( K a,l^aA)- (2.7) 

a 

For future purposes we also put a metric d : /too x Koo —■> [0,oo), and fix an origin 
0 G AToq. In proving our main results, the metric d does not play a major role. However 
for practical applications, such as determining the spectral dimension of the Schrodinger 
operator, one needs to understand the interplay between the metric d and the measure 
Poo', see Remark 2.9 and Section 6. 
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Let the potential V be a nonnegative, locally bounded measurable function on K 0O . 
(In general, V can be a real-valued, locally bounded measurable function which is bounded 
below. By adding a suitable constant to V one retrieves the case of a nonnegative poten¬ 
tial.) 

Assumption 2.5. There exists a self-adjoint Laplacian — A on L 2 (IL 00 , //q c ) [equivalently, 
a local regular Dirichlet form (£, F) on L 2 (A' oc , Too)]> and that the potential V(x) —> +oo 
as d{ 0, x) —> Too. 

Proposition 2.6. Under Assumption 2.5, the Schrodinger operator (—A T V), regarded 
as a sum of quadratic forms, is self-adjoint on L 2 (Koo, Too)? and has pure point spectrum. 

Proof. This uses the min-max principle as stated in [ 18, Theorem XIII.2], and then follows 
the proof of [48, Theorem XIII. 16]. □ 

By virtue of Proposition 2.6, we can define the eigenvalue counting function for (-A + 
V) on K^: 

N(Koo, Too, V, A) := # {A* (-A T V) < A} . (2.8) 

We are interested in the asymptotics of iV(iv 0O , Too, V, A) as A —> oo. In order to state the 
precise results, we will impose some mild conditions on the potential V. 

Given a potential V on K 0 Q , let V A (resp. V v ) be the function which is piecewise 
constant on each cell K a . and takes value sup xeA - fi V(x) (resp. inf xg x Q G(x)) on K a . We 
introduce the associated distribution functions 

F a (V,A) := ti 00 {{xeK 00 :V A (x)<\}), (2.9) 

F v (V,A) := Too ({z G Koo : V w (x) < A}) . (2.10) 

Note that F A (V, A) < F w {V, A). 

Assumption 2.7. There exists a constant C > 0 such that F V (G, 2A) < CF^{V, A) for 
all sufficiently large A. 

Note that this assumption implies that both A V (G, •) and T A (G, •) have the doubling 
property: there exist C v , C A > 0 such that 

F v (V,2A) < C w F w (V, X) and F A (V, 2A) < C A F A (G, A) (2.11) 

for all sufficiently large A. 

Assumption 2.8. The potential V on K 0 0 satishes 

F A (V] A) = 1 °(1) as ^ (2-12) 

Remark 2.9. To understand Assumption 2.7 or 2.8, it helps to keep the following example 

in mind. Let (K^, p 0 0 , d) be a metric measure space which admits a cellular decomposition 

into copies of the compact metric measure space (K,p,d). Let diam^(AT) be the diamater 
of K in the d-metric. Further suppose that Too is Ahlfors-regular: there exist positive 
constants c\, C 2 , and a such that 

cir Q < Too (B d (x,r)) < c 2 r a (2.13) 


7 



for all x E K oo and sufficiently large r > 0. As for the potential V, assume that there 
exist /3 > 1 and 7 E (0,1] such that 


C 3 d( 0 ,x)^ < V(x) < 04 ^( 0 , x)^, 


(2.14) 



(2.15) 


for all x,y E K^. By a direct calculation one can verify that (2.14) implies 


c 6 \ a/f} < F h (V, A) < c 7 A"/ /j , 


(2.16) 


which satisfies Assumption 2.7. Meanwhile, (2.15) implies 

E A (x) — E v (x) < c8[diamrf(A')] 7 d(0, x)^ -7 . 


(2.17) 


Thus (2.14) and (2.15) together imply Assumption 2.8. 


Our main results are the following. 

Theorem 2.10 (Existence of spectral dimension). Under Assumptions 2.1, 2.5, and 2.7, 
we have that 



(2.18) 


where F(V, A) := ({x E K : V(x) < A}). In particular, if F(V,X) = @(A^) as A — >■ 

00 , then d s (V) = + 2/3 is the effective spectral dimension of the Schrodinger operator 

(-A + V). 

Theorem 2.11 (Bohr’s formula). Under Assumptions 2.2, 2.5, and 2.8, 


lim NjKocPoo, V, A) 
A—^00 g(\ 5 , A) 


(2.19) 


where 


J [(A - V(x)) + ] ds/2 G Q log (A - E(x)) + ) 


ij(V, A) : 


Poo{dx), (2.20) 


and (/)+ = max{/, 0 }. 

In what follows we shall refer to g as “Bohr’s asymptotic function.” 

The proof of Theorem 2.11, discussed in Section 3, utilizes Dirichlet-Neumann brack¬ 
eting on the eigenvalue counting function and on Bohr’s asymptotic function. This is a 
relatively standard technique which is explained in the mathematical physics literature; 
see e.g. [ 18, §XIII], The novelty of our approach is to restate the sufficient condition on 
the potential V in terms of its distribution function, which allows us to extend the classical 
Bohr’s formula to a wider class of settings, such as on unbounded fractal spaces. 






2.2 Laplace transform version 

There are also analogs of Theorems 2.10 and 2.11 for the Laplace-Stieltjes transform of 
the eigenvalue counting function 

poo 

C(K 00 ,fi 00 ,V,t):=Tr Koo {e- t ^ A+v ^}= e~ xt N^, p^, V, dX). (2.21) 

Jo 

When V = 0 this is the trace of the heat semigroup associated with the bare Laplacian — A. 
More generally, it can be regarded as the trace of the Feynman-Kac semigroup associated 
to the Markov process driven by — A subject to kiling with rate V{x) at x £ Koo- 

The reason for stating the analog versions is because for certain compact metric mea¬ 
sure spaces, it is not known whether an explicit Weyl asymptotic formula for the bare 
Laplacian (Assumption 2.2) exists. However it may be the case that an asymptotic for¬ 
mula for the heat kernel trace (in some literature it is also called the partition function) 


C(K, fi, t) := Tr{e tA } = J pt(x,x) p(dx) (2.22) 

I< 

exists in the t f 0 limit. Here pt{x,y ) (t > 0, i,j/£ K) is the heat kernel associated to 
the Markov semigroup e t ^ generated by the self-adjoint Laplacian —A on L 2 (K, p). To 
be more precise, we denote by C h (K, p,t) the heat kernel trace of the Laplacian — A b on 
L 2 (K,p) with boundary condition b £ {A, V}. Then 

pOO p 

C h (K,p,t) = / e~ xt N h (K, p, dX) = / p\ fx, x) p{dx), (2.23) 

Jo Jk 

where N h (K, p, A) is as in (2.1), and pf(x,y ) is the heat kernel associated with the in¬ 
finitesimal generator — A b . 

Assumption 2.12 (Existence of the spectral dimension for the bare Laplacian). There 
exists a positive constant d s such that 

0 < lim t ds / 2 C b (K, fi,t) < lim t ds ^ 2 C h (K, y,t) < oo (2.24) 

40 HO 


for b £ {A, V}. 


Theorem 2.13. Under Assumptions 2.5, 2.7, and 2.12, we have that 


£(Koo, /Too; L, t) ..— T( U-oo , ^oo i L, t) 

< lim- -. 777- < 00, 

“ 40 t~ d °l 2 F(V,t) 


0 < lim . , 


(2.25) 


where 

T(V,t)= f e-^/ioo (dx). (2.26) 

Jl<oo 

In particular, if F(V, A) := ({2; £ Koo : V(x) < A}) = 0(A^) as A —>■ oo, then d s (V) = 

d s + 2/3 is the spectral dimension for the Schrodinger operator (—A + V). 


Assumption 2.14 (Weak Weyl asymptotics for the bare Laplacian). There exists a 
positive constant d s and a continuous function H : M + -A M+, independent of the boundary 
condition b £ {A, V} and with 0 < inf H < sup H < 00, such that as t j. 0, 


C b (K, p,,t) = t 


_ f—d s /2 


H(t) + p h (t) 


(2.27) 


where p h (t) denotes the remainder term of order o(l). 
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Theorem 2.15 (Laplace transform version of Bohr’s formula). Under Assumptions 2.5, 
2 . 14 , and 2.8, we have that 


l im C(K^, 

Moo ? V,t) 

m t- d */ 2 H(t)E{v, t ) 


(2.28) 


Note that (2.28) can also be interpreted as the asymptotic factorization of the trace 
of the Feynman-Kac semigroup: 


lim Tr Kx {e t( - A+y 
40 Tr A -{e iA } • Tr Koo {e~tV} 


= 1. 


(2.29) 


Remark 2.16. We make a few comments concerning the connections between Assumption 
2.1/2.2 and Assumption 2.12/2.14. 

(i) Assumption 2.1 is equivalent to Assumption 2.12. 

(ii) Assumption 2.2 implies Assumption 2.14. However, the reverse implication is possi¬ 
bly not true, since the classical technique of Tauberian theorems may not be appli¬ 
cable in this context. 

(iii) In order to prove Bohr’s formula (Theorem 2.11), we impose in Assumption 2.2 
that the function G be a periodic function. This is natural in light of the fractal 
examples we are interested in. However, to prove the Laplace transform version 
of Bohr’s formula (Theorem 2.15), one does not need to assume the log-periodicity 
in Assumption 2.14. This leads to the question of whether one could relax the 
periodicity of G and still be able to prove the original Bohr’s formula in greater 
generality (we do not address this question in the present work). 


2.3 Application of the main results 


To illustrate how our main results can be used, we now describe the “harmonic oscillator” 
problem on the Sierpinski gasket which was investigated in [ ]. For discussions of more 

general unbounded potentials on other fractal-like spaces, see Section 6. 


Example 2.17 (Harmonic oscillator on the infinite blow-up of the Sierpinski gasket). Let K 
be the Sierpinski gasket (SG). To construct SG, we first set the three vertices {pi,P 2 ,P 3 } 
of an equilateral triangle in M 2 , and then introduce the contraction maps : M 2 —y M 2 , 
^j{x) = \{x—Pj)+Pj, j = 1, 2, 3. Then SG is the unique fixed point K under the iterated 
function system consisting of the 'F 7 : K = Let w = w\W 2 • • • w m be a word of 

length \w\ = m where each letter Wj G {1, 2, 3}, and define the map o ■ ■ ■ o 

We endow SG with the uniform self-similar measure v with v(^> w K) = 3 - l w L The 
theory of Kigami [10] allows us to define the standard Laplacian on L 2 (SG, v) with either 
Dirichlet or Neumann condition on the boundary d(SG ) = {pi,P 2 ,P 3 }- Moreover, Kigami 
and Lapidus [33] proved that the eigenvalue counting function for the standard Laplacian 
satisfies 


JV b (SG>, A) 


yd a /2 



+ o{ 1 ) 


(b G {A, V}), 


(2.30) 


where d s = 2 log 3/log 5, and G is a cadlag periodic function with period ^log5 and 
contains discontinuities. Thus Assumption 2.2 is satisfied. 

Next, for each infinite word w = W 1 W 2 ■ ■ ■ which is not eventually constant, define 


SGZ:= U Ko---o^)(5G) (2.31) 

m=0 
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to be the infinite blow-up of SG associated with the word w. This is an unbounded fractal 
space where the neighborhood of any point x € Koo is homeomorphic to SG, and thus 
is a fractal analog of a manifold, called a fractafold by Strichartz [56]. Properties of the 
Laplacian on SG^fi are discussed in [ ]. Here we point out that by construction, SG ™ 

admits a cellular decomposition into copies of SG which intersect on the boundary only. 
Thus the measure v on SG can be readily extended to the measure Vqq on SG ™ . 

In [19] Fan, Khandker, and Strichartz studied the spectral problem of a harmonic 
oscillator potential V on a class of infinite blow-ups of SG. They defined V to be a 
solution to — AV = —1 on SG “ which grows unboundedly as d(0,x) oo and attains 
a minimum at some vertex xo £ Koo- (The first condition is a suitable replacement of 
V(x) = |x| 2 , which is available only in the Euclidean setting.) Note that this implies that 

V(x) grows at infinity at rate comparable to a positive power of R(xo,x), where R(-,-) 
is the effective resistance metric on SG^fi. This verifies Assumption 2.7. However we 
cannot verify Assumption 2.8 for general words w. Paper [L9] also contains information 
about spectral dimension, which depends on the blow-ups of SG. Through a mix of 
computations and numerical simulations, the authors of [19] were able to find properties 
of the low-lying eigenfunctions, as well as the asymptotic growth rate of the eigenvalue 
counting function of —A + V [19, Theorem 8-1 and Eq. (8.18)]: 

cX d ° < N{SGZ ,/Zoo,E,A) < C\ da . (2.32) 

Among the open questions posed in [19, Problem 8-3 and Conjecture 8-4] is finding the 
asymptotic “Weyl ratio” A _ds ^^ 2 A'(lF 0O , /i,-*,, V, A) of the eigenvalue counting function. 
Here we can provide an indirect answer. Given that Assumptions 2.2, 2.5, and 2.8 are 
satisfied, Bohr’s formula (Theorem 2.11) says that as A —> oo, 

N(SG% 0 ,v oo ,V,X) = (l + o(l)) / [(A-E(x)) + ] dfl/2 Gf^log(A-H(x)) + ') d^(x). 

(2.33) 

This in some sense answers the Weyl ratio question, in spite of the non-explicit nature of 
the integral on the right-hand side. 

The rest of this paper is organized as follows. In Section 3 we describe the tools needed 
to establish Bohr’s formula in the setting of an unbounded space which admits a cellular 
decomposition according to the setup in Section 2.1. In Section 4 we show how to restate 
the general sufficient condition for Bohr’s formula in terms of distribution functions of V v 
and V A , and also give a “weak” version of Bohr’s formula. We can show how the addition 
of an unbounded potential leads to the absence of gaps in the spectrum of — A + V. This 
is of independent interest since the spectrum of the bare Laplacian on certain fractals ( e.g . 
the Sierpinski gasket) has gaps. In Section 5 we establish the Laplace transform version 
of Bohr’s formula. Finally, in Section 6, we discuss applications of our main results to 
various unbounded potentials on several types of unbounded fractal spaces. 

3 The general Bohr’s formula 

In this section and the next section, Assumptions 2.2 and 2.5 are in force. 

3.1 Bohr’s asymptotic functions 

Let — A A (resp. — A v ) be the Laplacian on ifi^Koo, fioc) with Dirichlet (resp. Neumann) 
conditions on the gluing boundary U a dK a . For each potential V, let V A (resp. V v ) be 
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the piecewise constant function which takes value sup xeft - ri V(x) (resp. inf x &K a V(x)) on 
K a . Thanks to Proposition 2.6, one can introduce the eigenvalue counting functions 


NiKocfi oo,^, A) := #{A*(-A + y)< A}, (3.1) 

■= #{X i (-A A + V A ) < A}, (3.2) 

^(iWocu^A) := # {A* (-A v + V w ) < A} . (3.3) 

Note that since (—A v + V w ) < (—A + V) < (—A A + V A ) in the sense of quadratic forms, 

N A (K 00 ,fi 00 ,V,X)<N{K 00 , f i 00 ,V,X)<N v (K 00 , f i 00 ,V,X). (3.4) 


We shall show that under some mild additional conditions on V, iV(A' oc , V, A) is 
asymptotically comparable to the “Bohr’s asymptotic function” 

g(V, A):= j [(X-V(x)) + ] ds/2 G(^log(X-V(x)) + ^ dti oo(x), (3.5) 


where (/)+ := max{/, 0}, and G is as appeared in Assumption 2.2. In order to estimate 
this rate of convergence, we introduce the functions 





(3.6) 


and 


n h (V, A) := f [(A - V h (x))+] ds/2 R h ((A - F b (x))+) d^x) 

K oo 


(3.7) 


for b 6 {A,V}, where R h is the remainder term which appeared in Assumption 2.2. 
Observe that since V h (x) is constant on cells, the right-hand side expressions in (3.6) and 
(3.7) are really discrete sums: 


9°(V, A) 
77 b (V,A) 


E 

{a:V b \ <A} 

1 'Iq 

E 

{a:V b \ <A} 

1 IXQ- 


A - V* 


X — V t 


K a 


K a 


d a /2 

Gl- log(A-W 


K a 


ds/2 


R b A - V t 


K c 


(3.8) 

(3.9) 


Moreover, by Proposition 2.4, decouples into the various K a according to the Dirichlet 
or Neumann boundary condition, so 


A rb (/i 00 ,^ 00 ,P,A) 


E " b 

{or-V* | <A} 

I i\.Ql 


(k<x, 9a, X-V h 



Pulling (2.3), (3.8), (3.9), and (3.10) together we obtain 

N h (K 00 , t i 00 ,V,\) = g b (V,\)+K b (V,\). 


(3.10) 


(3.11) 
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3.2 Monotonicity of Bohr’s asymptotic functions 

A key monotonicity result we need is 

Proposition 3.1. Fix a potential V. Then each of the functions A i —> g(V, A), A i —> 
g A {V, A), and A g w (V, A) is monotone nondecreasing for all A > 0. Moreover, g A (V, A) < 
g{V,X)<g\V,X). 

This follows from the monotonicity of the integrand of the g function. 

Proposition 3.2. The function W(X) = X ds ^ 2 G log A) is nondecreasing. 

Remark 3.3. This proposition implies, in particular, that G has a cadlag version. Although 
the result is very simple, we could not find it in the literature. 

Proof of Proposition 3.2. If W was not nondecreasing, then there would existed A2 > Ai > 
0 such that W(X 2 ) = VF(Ai) = — 5 < 0, which contradicts (2.3) and the monotonicity of 
N h {K,g,X). O 

Proof of Proposition 3.1. Fix a potential V. For each A > 0 and x £ K^, put 

IF(A, V, x) = ((A - V(x)) + ) d ‘/ 2 G Q log((A - U(x)) + )) (3.12) 

and 

W h (X, V. , x) = ((A - V h {x)) + ) d ‘/ 2 G Q log((A - P b (*)) + )) • (3.13) 

Observe that W(X,V,x) = W({X - V(x))+) and W h {X,V,x) = W{(X - P b (x))+). Us¬ 
ing Proposition 3.2 we deduce the following two consequences. First, A e-)- W(X,V,x) is 
nonnegative and monotone nondecreasing for each x. And since g(V, A) is the weighted 
integral of W (A, V, x) over x, it follows that A 1 —> g(V, A) is also monotone nondecreasing. 
The monotonicity of A 1— > g h (V , A) is proved in exactly the same way. Second, the mono¬ 
tonicity of W( A) implies that VF A ( X,V,x) < W(X,V, x) < W w (X,V, x) for each x, and 
upon integration over x we get g A (V, A) < g(V, A) < g v (V, A). □ 

3.3 Bohr’s asymptotics via Dirichlet-Neumann bracketing 

We have all the necessary pieces to state the error of approximating N(K^, V, A) by 

g(v, A). 

Theorem 3.4 (Error estimate in Bohr’s approximation). Under Assumptions 2.2 and 
2.5, we have 


N(I <ocMocb A) 


- 1 


< max 
b£{A,V} 


g h (V,X) _ i + n h (V,X) 


g h (v, A) 


5 b (^,A) 


g(v, A) 

where b = A (resp. b = Vj if b = V (resp. if b = A). 

Proof. From (3.4) we have 

V a (A^oo, Moo, V,X)<N(K 00 ,g 00 ,V,X)< V v (A to , Moo, V, A). 


(3.14) 


(3.15) 
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Meanwhile by Proposition 3.1 


g A (V,X)<g(V,X)<g^V,X). 


(3.16) 


Therefore 

N A (K 00l lioo, V, A) ^ IV (Koodoo, V, A) ^ /*«*,, V, A) 

g v (V, A) - g(V, A) - ^(v.A) • 1 • j 

Subtract 1 from every term in the inequality (3.17), and then use (3.11) to replace 
A rb (iV 0O , /Tqo, V, A) with g h (V, A) + 77 6 (V, A). Finally, we can estimate the absolute value 
of the middle term of the inequality by the maximum of the absolute value on either side 
of the inequality. □ 


Having established the main error estimate, Theorem 3.4, we can now give an abstract 
condition on V for which Bohr’s formula holds. 

Assumption 3.5. The potential V on Ah*, satisfies 

= 1 + o(l) as A —> oo. (3.18) 

Theorem 3.6 (Strong Bohr’s formula). Under Assumptions 2.2, 2.5, and 3.5, we have 


lim 

A—>oc 


N(K 00,/Xoo, V, A) 
g(V, A) 


= 1. 


(3.19) 


Proof of Theorem 3.6. Observe that Assumptions 2.2 and 3.5 together imply that the 
error term stated in Theorem 3.4 is o(l). □ 


4 Connection between Bohr’s formula and the 
distribution function of the potential 

Assumption 3.5 can be too abstract for applications dealing with fractal spaces. We now 
explain how this assumption can be restated in terms of distribution functions of V: 

F(V, A) := ^({x € A'^ : V(x) < A}) and F h (V, A) := ^{{x € A TO : V h (x) < A}. 

(4.1) 

Lemma 4.1. We have that 

rW( A) rW(X) 

g(V, X)= F(V,X-W~\t))dt and g h (V,X)= F h {V,X-W~ 1 (t))dt, 

Jo Jo 

(4.2) 

where 

VF _1 (t) = inf {A > 0 : W{X) > t} (4.3) 

is the generalized inverse ofW(X) = X ds ^ 2 G(^ log A). 


14 







Proof. We start with a fundamental identity in measure theory. For any nonnegative 
function / on a a -finite measure space (X, m), Fubini’s theorem tells us that 



f(x) m(dx) 



m({x G X : f(x) > f}) dt. 


Applying this identity to g(V, A) we find 


(4.4) 


r r oo 

9 (V,X)= / W((X-V(x))+)d f i 00 {x)= / gocdxeK^: W((\ - V(x))+) > t}) dt. 

J Koo J 0 

(4.5) 

Since W is monotone nondecreasing (Proposition 3.2), it has a well-defined generalized 
inverse W ~ l , which satisfies 


W(X) > t <*=>■ A > IF _1 (f). (4.6) 

So the right-hand term in (4.5) can be further rewritten as 

/Uoo ({* € Koo : (A - V{x))+ > W-\t)}) dt. (4.7) 

Now by assumption V is a nonnegative potential, so W~ l {t) < (A — V(x))+ < A, or 
equivalently, t < W(A). This places an upper bound on the integral, and we get 

rW{ A) rW( A) 

g(v, A) = / Hoo ({x g Xoo : V{x) < A - dt = / F(V, A - W-\t)) dt. 

Jo Jo 

(4.S) 

The proof for g h (V, A) is identical. □ 



Observe that for X < X 1 , 

g w (V,X)-g A (V,X') = 


and 


r-W(A) 


[F v (V, X - W~ l {t )) - F a (V, A' - W~\t))\ dt 


rW(X') 

/ F A {V,X' -W-\t))dt, 
Jw( A) 


(4.9) 


fW( A) 


g v (V,X')-g\V,X) = / [-F v (V, A' - W _1 (t)) — P A (V, A — JT _1 (t))] dt 


+ 


rW( A') 
lw( A) 


F v (V,A' - W” 1 ^))^. 


(4.10) 


These identities suggest that if the difference of the distribution functions F v (V, A) — 
F A (V, A) can be controlled, then one can control the difference g v (V, A) — <? A (V, A). Indeed 
we have 


Proposition 4.2. Assumption 2.8 implies Assumption 3.5. Therefore, the strong Bohr’s 
formula (Theorem 3.6) holds under Assumptions 2.2, 2.5, and 2.8. 
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Proof. Let h(V, A) = - 1 > 0 

0 < <7 v (V,A)-</ a (V,A) 

rW( A) 

= / [l + h(V,X-W~\t)) -1}F A (V,X-W~\t))dt 


). Then 


< 


( sup h(V,X-W-\t))\ r w 
\o<t<w(X) ) Jo 


F A (V,X-W~\t))dt 


sup h(V,s) g A {V, A). 
, 0<s<A / 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


Assumption 2.8 implies that sup 0<s<A h(V, s) = o(l) as A —> oo, so we obtain Assumption 
3.5. “ “ □ 


4.1 A weak version of Bohr’s formula 


Motivated by [19,44,45], we also give a weak version of Bohr’s formula as follows. 
Theorem 4.3 (Weak Bohr’s formula). Let X* > X with X* — A = o(A) and 


F v (V, A) 


= 1 + o(l) and 


F A {V, A) 


F A (V, A*) v ' F V (R,A* 

Then, with Assumptions 2.2 and 2.5, we have 


= l + o(l). 


(4.15) 


lim 

A^-oo 


N^Kqq, poo, V, A) 

ff(v, X*) 


(4.16) 


The statement of Theorem 4.3 is reminiscent of the situation when one compares two 
nondecreasing distribution jump functions with closely spaced jumps. When the jumps 
asymptotically coincide, then the difference of corresponding measures tends to zero in 
the sense of weak convergence. 


Proof. By mimicking the proof of Theorem 3.4 we get 


N(Koo,fj,oo, V, A) 

g(v, x*) 


< max 
b£{A,V} 


g S (KA) 

g h (V,X*) 


1 + 


ft B (V, A) 
9 h (V,X*) 


(4.17) 


Since A* — A = o(A) as A — > oo, the ratio 77 b (V, A )/g h (V, A*) can be made to be o(l). So 
the key estimate is to show that g h (V, X)/g h (V, A*) = 1 + o(l) for both b = A and b = V. 
(This is to contrast with the case X' = X as shown in Proposition 4.2, where a one-sided 
bound suffices because g w (V, A) — g A {V, A) > 0.) 

From (4.9) we find 


\g w (V,X)-g A (V,X*)\ < W(X) ( sup [F^(V, s) - F A {V, s + A* - A)]] ) (4.18) 

\0<s<\ / 

+[W(A*) — W(X)\ ( sup F a (F,s)| (4.19) 

yO<s<A*—A J 

According to the first condition in (4.15), sup 0<s<A [F v (F, s) — F A (V,s + A* — A)]] = 
o(F w (V, A)) and sup 0<s<A *_ A F A (V, s) = o(F A (V, A*)). This implies that the absolute 
value on the RHS of (3.14) is o(l) for b = A. Similarly, the second condition in (4.15) 
implies that the absolute value on the RHS of (3.14) is o(l) for b = V also. □ 
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5 Laplace transform (heat kernel trace) version 
of Bohr’s formula 

In this section we impose Assumption 2.5 and either one of Assumptions 2.12 and 2.14, 
and prove Theorems 2.13 and 2.15. Let us introduce the traces 


£(1Loo,Moo, V,t) 
{Koodoo, V,t) 
C w (ifoc/ioo, V,t) 


= Tr Koo {e-^- A+v ^}, 

= 'lWe-* ( - AA+yA) } 

-t(-A v +y v 


= Tr ^oc{e" 




Observe that C A {Koo, /x^, V, t) < £(iLoo, [Xoo, V, t) < £ v (, fXoo , V, t). 
Since L 2 (A' 00 , /Xoo) = 0 a L 2 (K a , /x a ), it follows that 

£^# 00 ,^ 00 , V,t) = ^ C h (K a ,fx a , V, t), b € {A, V}, 


where 


C 


\K a , n a , V, t) = Ti Ka (e t( Ab + yb )| = C h {K a ,Ha,t) - exp (^-t U b ^ ^ 


Let 


(5.1) 

(5.2) 

(5.3) 


(5.4) 

(5.5) 

(5.6) 


(5.7) 


T(V,t) := / e~ tv W Hoc (dx). 

Similarly define 

T h (V,t):= f e- tv "W fi^idx) 

for b 6 {A,V}. Observe that J rA (V,t) < F(V,t) < J rV (V,t), and that Assumption 2.5 
ensures that iFfV, t) and J rb (V, t) are finite for t > 0. 

Proof of Theorem 2.1 3. Let us first note that 


By (5.4), 


C A {K00,1X00, V,t) < C{Kqo,IXqo, V,t) < C V (Koodoo, V,t) 
t -ds/^^(V, t ) - t-^/ 2 A(F,t) “ t~ d ^ 2 F A {V,f) 


C h {Ko o, /Xoc, V,t) = J 2 £ h {Kon fia, V., t) 

a 

= ^2 £ b (K a ,haj t ) • exp ( -t U 1 
a 

= ^£ h (K a ,fx a ,t)- f 


K 0 


a -tV h (x) 


/X a (dx) 


(5.8) 

(5.9) 

(5.10) 

(5.11) 


K a 


Under Assumption 2.12, there exist positive constants C\ and C 2 such that for all suffi¬ 
ciently small t, 

t d ‘/ 2 Py (Koodoo, V,t) < OiV/ e~ tvW ^ f x a (dx) = C 1 ^(V,t), (5.12) 

a J K a 

t ds/ 2 C A (Koodoo, v,t) > C 2 J2[. e~ tvA W fx a (dx) = C 2 F A (V,t). (5.13) 
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Meanwhile, by Fubini’s theorem and by the nonnegativity of V, we have 


' — OO 

r oo 


4oo ^ 

{ X G Koc 

. e -« /l> p) > 

(5.14) 

4oo ( 

G Koc 

, : e“ tvb ( x) > e" a ) te~ tx dX 

(5.15) 

Too ^ 

G i^oo 

: V h (x) < A}^ te~ tx dX 

(5.16) 

F b (V, X)te~ tx 

dX. 

(5.17) 


Hence under Assumption 2.7, there exists Aq > 0 such that 


F v (V,t) = I F w {V, X)te~ tx dX 


(5.18) 

(5.19) 

< F v (V, A 0 ) J* te~ tx dX + C j F A (v, te~ tx dX (5.20) 

'(V,A 0 ) (l-e~ tXo ) +C f F A {V 1 X)-2te~ 2tx dX (5.21) 

V 1 j An/2 


1 0 

f\o roo 

/ F v (V, X)te~ tx dX + / F v (V, X)te~ tx dX 
1 0 J A 0 

f Ao 


= F x 


Therefore 


< F v (V, A 0 ) (l - e" a °) + C.F a (V, 2 1). 




F A (V,t) F A (V,2t) " ~ ' * K ' , '" J/ F A (V,2t)' 

Since lirri^o(1 — e _tA °) = 0 and t1->- F A (V, 21) is monotone decreasing, it follows that 


(5.22) 

(5.23) 


40 F A (V,t ) 
Putting everything together we find 


ttv (y 4-\ _1 _ p -t\ 0 

lim X <C + F v (V, A 0 ) lim = C. 


40 F A (V,2t) 


— £(^oo,/ioo, V,t) ( p— t ds / 2 £ V (Kqq , ^oo , V, f) \ / ^ V (P,f) 

/o-r-ZT, -s - I ™ JV4oJ zrA (P,f) 


40 t~ d °/ 2 F(V,t) ~ y 40 A v (V,t) 

.. CiK^fiooM) ^ 0 ,V,i)\ A A (V,t) 


lim 


— t -d s /2jr {Vit) 


> lim 


40 ^ A (V,i) 


lim 


40 ^G 7 ,*) 


Thus 


CaC- 1 < lim £ GWoo,*) < HE *2 < CiC. 


(5.24) 

(5.25) 

(5.26) 

(5.27) 


40 t~ d °/ 2 F{V,t) ~ 40 t~ d ‘/ 2 F(V,t) 

Finally, regarding the spectral dimension of — A+H, we note that F(V, A) = 0(A' 3 )a-> O o 
is equivalent to F(V,t) = 0(i~ ,3 )4o> an easy consequence of Laplace transform. Thus 
according to (5.27), £(Koo, ^oo, V,t) x £-( rf s+2/3)/2 as £ ^ q. □ 


Proof of Theorem 2.15. Combining (5.11) with Assumption 2.14 we obtain 
^(iLoo^oo, V,t) = |ff(i) + ^(i)! A b (V,t) 


(5.28) 
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which, together with (5.8) and after some manipulation, implies 


CjKoo, Moo? v,t) 
t- d '/ 2 H(t)F(V,t) 


< max 
b£{A,V} 


(, | pHt)\ 
[ + #(t) J 


(5.29) 


Next, by Assumption 2.8 and (5.18), 


■F v (V,t) 


roc 

F A (V, t) + / (F V (H, A) - F A (H, A)) te~ tx d\ 

Jo 

F A (V,t) + o(F A (V,t)) m 


(5.30) 

(5.31) 


because for any 5 > 0 there is A<5 > 0 such that F V (F, A) — F A (V, A) 
A > A$. Thus = 1 + o(l) as 1 } 0. Hence (5.29) implies (2.28). 


< <5F v (V,A) when 
□ 


6 Examples 

In this section we provide several instances on both classical and fractal settings whereby 
the existence of the spectral dimension of —A + V can be proved, and moreover, Bohr’s 
formula holds. 


6.1 Euclidean spaces 

One would be remiss not to mention the most classical setting, which is the Schrodinger 
operator — A + V on M rf , where A = Yli=ii^ 2 / ^ x l) an d V is an unbounded potential. See 
e.g. [48, §XIII. 15]. The key idea is to partition (the unbounded space K^) into cubes 
of side 1 (the cells K a ). Then by applying the machinery outlined in the previous section, 
one arrives at the following well-known result: if V(x) = 0(|x|^) as |x| oo, then Bohr’s 
formula holds, and the spectral dimension of this Schrodinger operator is d( 1 + 2//3). 

In dimension 1 Bohr’s formula can be established for logarithmically diverging poten¬ 
tials. The proof method involves solving a Sturm-Liouville ODE, which appears rather 
particular to one-dimensional settings, and may be difficult to generalize to higher dimen¬ 
sions. We refer the reader to [28,37] for more details. 

6.2 Infinite fractafolds based on nested fractals 

Nested fractals are introduced first by Lindstrpm [ ]. The typical examples to keep in 

mind are the Sierpinski gaskets SG(n), where n denotes the length scale of the subdivision. 
There are also higher-dimensional analogs of SG. 

On nested fractals, and more generally post-critically finite (p.c.f.) fractals, one can 
define a notion of the Laplacian (or a Brownian motion). See e.g. [5, §2~§4], [40, Chapters 
2~3], [58, Chapters 1~2] for the relevant definitions and results. We will need just one 
result on the spectral asymptotics of the Laplacian on p.c.f. fractals with regular harmonic 
structure. 

Proposition 6.1 ([33, Theorem 2.4], [40, Theorem 4.1.5]). Let K be a p.c.f. fractal, 
and p, be a self-similar measure on K with weight {ni)f =1 - Assume that pnri < 1 for all 
i E (1,2, , IV}. Let d s be the unique number d which satisfies YliLilt = C where 

7 i = yjrifii. Let N A (K,fi, A) (resp. N w (K,fi, A)J be the eigenvalue counting function for 
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the Laplacian on L 2 (I\,p,) with Dirichlet (resp. Neumann) boundary condition. Then for 
b € {A, V}, 

0 < lim X~ ds/2 N h (K,ii,X) < BE X~ ds/2 N h (K, p, X) < oo. (6.1) 

A—>oo A ^°° 


Moreover: 

(a) Non-lattice case: If YliL i ^ l°g hi a dense subgroup o/M, then the limit 

lim X~ ds/2 N h (K,n,X) 

A —f oo 


exists, and is independent of the boundary conditions. 

(b) Lattice case: ^ log 7i * s a discrete subgroup o/M, let T > 0 be its generator. 

Then as A —>■ oo, 


N h (I\, n, X) 



yds/ 2 


( 6 . 2 ) 


where G is a right-continuous, T-periodic function with 0 < inf G < sup G < oo, 
and is independent of the boundary conditions. 


We remark that the proof of Proposition 6.1 relies upon Feller’s renewal theorem [20]. 

Our goal is to state Bohr’s formula for the Schrodinger operator on a class of unbounded 
fractal spaces. As mentioned in Section 1, one such candidate is a fractafold based on a 
nested fractal. We shall consider two types: 

(i) The infinite blow-ups of a nested fractal in M d , d> 2. (See Figure 1). 

(ii) Infinite periodic fractafolds based on the planar Sierpinski gasket K = SG(n), 
equipped with a metric R. (In practice, R is taken to be the resistance metric, but 
the results to follow do not depend explicitly on the specifics of R.) The examples 
we will consider can be constructed by first defining an infinite “cell graph” T, and 
then replacing each vertex of T by a copy of K, and gluing the K a in a consistent 
way. With this construction the metric R on K extends to a metric R on K ^ in the 
obvious way. For instance, one can construct the ladder periodic fractafold (Figure 
2) and the hexagonal periodic fractafold (Figure 3). 

To establish Bohr’s formula, we will need information about the measure growth of 
balls in K^. For the infinite blow-ups of a nested fractal, it is direct to verify that for all 
x E Koo and r > 0, 

cr df ’ R < Hoo(Bfi(x,r)) < Cr df ’ R , (6.3) 

where df t R is the Hausdorff dimension of K with respect to the metric R on K. 

For the periodic fractafolds a slightly different analysis is needed. Let dr be the graph 
metric of the cell graph T, and Bd r (z,r ) := {y E T : da(z,y) < r} be the ball of radius 
r centered at z in T. Since I \is constructed by replacing each vertex of T by a copy of 
K , we can estimate the volume growth of balls in using the cardinality of balls in F. 

Proposition 6.2. Let D(I \) := diamR(AT). For all x E K ^ and all r > 2 D{K), 


Too ( B dr r - 2 D(K))) < Too{Br(x, r)) < ( B dv ('if(x), r + 2D{K))) , (6.4) 

where if(x) is the vertex in T which is replaced by the cell K a 3 x in the periodic fractafold 
construction. 
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Proof. Let ?/(r) := r/D(K) > 2. Then Br(x,v) = Bjt(x,r](r)D(K)) and 


B R (y,([ V (r)\-l)D(K))CB R (x,r,(r)D(K))CB R (y,(\r,(r)] + l)D(K)) (6.5) 

for any y which lies in the same cell K a as x. Here |_aj (resp. [a]) denotes the largest 
integer less than or equal to a (resp. the smallest integer greater than or equal to a). It 
is then direct to show that there exist y such that Br(p, (["77(7-)] T l)D(K)) is covered by 
the union of all cells I\ a which are at most distance (|b?(r)] + 1) from y in the T metric. 
Since each cell has /x-measure 1, the /x-measure of the cover is equal to the cardinality of 
Bd T {^{x), [^(r)] T 1). The upper bound in (6.4) follows by overestimating |b/(r)~| + 1 by 
r/(r) + 2. The proof of the lower bound is similar. □ 

We can now state the main result of this subsection. 

Proposition 6.3. On the infinite blow-up of a nested fractal (resp. the ladder periodic 
fractafold based on SG(n), the hexagonal periodic fractafold based on SG(n)), Bohr’s 
formula holds for potential of the form V{x) ~ R(t),xf for any (3 > 0. In particular, the 
spectral dimension of —A + V is d s (V) = d s + 2 {dh/(3), where dh equals the Hausdorff 
dimension of the nested fractal with respect to the metric R (resp. 1, 2). 

Proof. Since each K a which makes up the cellular decomposition of is isometric to 
the same nested fractal I\, by Proposition 6.1 we have that Assumption 2.2 holds. 

Because the cells I\ a intersect at boundary points in a natural way, the Dirichlet form 
{£,3F) corresponding to the Laplacian —A on L 2 (IL 00 , /Xqo) can be built up as a sum of 
the constituent Dirichlet forms on L 2 {K a , p a ). Hence one can show self-adjointness of 
—A in the sense of quadratic forms. And since the potential V(x) grows unboundedly as 
d(0,x) —> Too, Assumption 2.5 then implies that (—A T V) has pure point spectrum. 

For condition (i), one can confirm that there exist constants c and C such that for all 
x 6 AToo and all sufficiently large r > 0, 

cr dh < Hoc[B r{x,v)) < Cr dh . (6.6) 

For the infinite blow-up (6.6) follows from (6.3) with d = dh r. As for the periodic 
fractafolds, note that the corresponding cell graphs T satisfy 

\Br(z, r)| x r dh ’ T for all z G T and r > 0, (6.7) 

where dh.r equals 1 (resp. 2) in the case of the ladder fractafold (resp. the hexagonal 
fractafold). Combining this with Proposition 6.2 we get (6.6) with dh = dh r- In all cases, 
we the find 


F( A) = /Xoo({* : V(x) < A}) ~ /ToCMO, A 1 ^)) ~ \ dh ^, (6.8) 

and the same asymptotics holds for F A (A) and F V (A). Finally, to see that condition (ii) 
holds, we use the inequality 

[F a (x)-F v (x)] < [R(0,x) + lf-[R(0,x)~ if <Cg[R{0,x) + if' 1 , (6.9) 

where Cp is an explicit constant depending on (j only. Observe that the RHS is uniformly 
bounded from above by a constant multiple of A 1_/3 1 for all x in the set {x : R v (x) < 
A}. □ 
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Figure 5: The double-ladder fractal field based on SG( 2). 

6.3 Infinite fractal fields based on nested fractals 

There is another notion of an unbounded space based on compact fractals, which are 
known as fractal fields. The name originated from Hambly and Kumagai [25], who 
were interested in studying fractal penetrating Brownian motions. Fractal fields differ 
from the fractafolds of the previous subsection in that we do not require neighborhoods 
of (junction) points in K^ to be homeomorphic to K. 

First consider the triangular lattice finitely ramified Sierpinski fractal field introduced 
in [60, §6], see Figure 4. Notice that this fractal field admits a cellular decomposition 
whereby cells adjoin at boundary vertices of SG(n). As a result, the proof strategy from 
the previous Proposition 6.3 applies in this setting. 

Proposition 6.4. On the triangular lattice finitely ramified fractal field based on SG(n), 
Bohr’s formula holds for potential of the form V(x) R( 0, xf for any f3 > 0. In partic¬ 
ular, the spectral dimension of —A + V is d s (V) = d s + (4 //3). 

Next we consider the double-ladder fractal field based on SG( 2), see Figure 5. An 
important difference here is that pairs of SG( 2) cells may adjoin either at a point or 
along a boundary segment, which makes this space infinitely ramified. In order to analyze 
this example using our methods, one needs to understand the eigenvalue problem for the 
Laplacian on SG(2)q := SG( 2) \ dl 2, where dQ consists of the top vertex and the bottom 
edge of SG( 2). This was investigated by Qiu [46], whose result we quote below. 

Proposition 6.5 ([46, Theorem 3.10]). Let IVq(A) be the eigenvalue counting function 
for the Laplacian on SG( 2 )q with boundary condition b G {A, V} on the top vertex and 
the bottom edge of SG(2). Then there exists a cadlag log 5 -periodic function G : M — > M, 
with 0 < inf G < sup G < oo and independent of b, such that 

Nq(X) = G(logA)A log3 / log5 + 0 (A log2 / log5 logA) (6.10) 

as A —> oo. 

Using this result we can show the validity of Bohr’s formula in this setting. 

Proposition 6.6. On the double-ladder fractal field based on SG{2), Bohr’s formula holds 
for potential of the form V (x) R{0,x)P for any fi > 0. 

Proof. The one nontrivial assumption to check is Assumption 2.2, which is furnished by 
Proposition 6.5. The other two assumptions, 2.5 and 2.8, are verified easily. The result 
then follows from Theorem 2.11. □ 
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Figure 6: The SG( 2) fractal field (or gasket tiling) considered in [25]. 


There are some obvious extensions of the double-ladder fractal field example, which 
we leave to the reader. An interesting open problem is to study the applicability of Bohr’s 
formula to the original fractal field (or gasket tiling) in [ ], shown in Figure 6. We note 

that heat kernel estimates are established on this fractal field [25, Theorem 1.1]. However, 
to the best of the authors’ knowledge, there is no corresponding Weyl asymptotic (or heat 
kernel trace asymptotic) estimate which is sharp to an o(l) remainder. In particular, the 
fact that the SG cells adjoin along edges rather than at points makes the analysis more 
delicate. 

6.4 Infinite Sierpinski carpets 

Let F C (d > 2) be a generalized Sierpinski carpet in the sense of [6,7], and let F n 
be its nth-level approximation. Following [6], we call F = |J ngN the pre-carpet, and 

Foe = UneN 0 the infinite carpet. The difference between the two is that F is tiled 
by unit squares and has nonzero Lebesgue measure, whereas F^ is tiled by copies of the 
same Sierpinski carpet F and has zero Lebesgue measure. In both cases, we adopt the 
Euclidean metric | • | and regard (K 0 c ,Foo, | • |) as the metric measure space, which has 
volume growth 

c\r df < Hoo(B(x,r)) < c 2 r df (x € A'oo, r > 0), (6-11) 

where df = (logm/logf) is the Hausdorff dimension of the carpet F with respect to the 
Euclidean metric. 

Proposition 6.7. Bohr’s formula holds on the pre-carpet F with potential V(x) ~ |x|^ 
for any fi > 0. In particular, the spectral dimension of (—A + V) on F is d + 2 (df//3), 
where d is the dimension of the ambient space M. d in which F lies. 

The case of the infinite carpet is more nuanced. Hambly [24] and Kajino [38] proved 
that the heat kernel trace of the bare Laplacian on F satisfies Assumption 2.14, with H a 
continuous periodic function of log t (though it is NOT known whether H is non-constant). 
Kajino [39] further showed the asymptotics of the heat kernel trace to all orders of the 
boundary terms. Note that their results imply that the eigenvalue counting function 
satisfies the asymptotics ciA rfs / 2 < N h (F, fi, X) < C 2 \ ds ^ 2 , but do NOT necessarily imply 
the sharper estimate, Assumption 2.2. As mentioned earlier, this is because the classical 
techniques of Tauberian theorems cannot be applied here. 
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Proposition 6.8. The Laplace transform version of Bohr’s formula holds on the infinite 
carpet F 00 with potential V{x) ~ \x\^ for any (3 > 0. In particular, the spectral dimension 
of (-A+V) on Foo is d s +2(df / /3), where d s is the spectral dimension of the bare Laplacian 
on F. 

Proof. By [24, Theorem 1.1] and [38, Theorem 1.2], Assumption 2.14 is satisfied on the 
constituent Sierpinski carpet F. In fact, [39, Theorem 4.10] provides a sharper result of 
the form 


a 

£ b (F,/i,t) = t~ d ^ 2 H{-logt) + ^t- dk / d -G\{-\ogt) + 0 (expt-cfT^r)) (6.12) 

fc=l 

as t f 0, where H and the G\ are continuous periodic functions, dk is the Minkowski 
dimension of F n {x = (x\,--- ,Xd) 6 : x\ = • • • = Xk = 0}, and d s and d w are 

respectively the spectral dimension and the walk dimension of F. 

We turn our attention next to the potential term F h (V,t). It is direct to verify that 
for any /3 > 0, 


J e dpoo^x) < J 


°° e -t\ dhoc({x : \x\P < A}) dX 


dX 


(6.13) 


K a 


poo 

= / te- tx Hoo{B{Q,\ 1,P ))d\ (6.14) 

Jo 

poo 

< c 2 t e~ tx \ df /P d\ < C 2 (df, /3)t df /P, (6.15) 

Jo 


and similarly 


e~^ x ^ diioc(x) > C 1 (d f ,/3)t df ^. 


(6.16) 


K a 


Using the inequality e s > 1 + s for s G M, we find 


e -t\x-y\P _ e -t\x-z\P 


< max I 


(e~t\x-y\\ e -t\x-z\^ . t ^ x _ y \h _ \ x _ z fj 
< Cp-t- max (e-M*-v\\ e -t\x-z\ p 'j . | y _ z |. 


It follows that as t f 0, 

F w (V,t) - F A (V,t) < C ■ 0(tF(V,t)) = o(F(V, t)), 
leading to the error estimate 

£(-F°o,Moo, V,t) 


- 1 


_ q do-di)/d w ^ 


t~ d sFH(-\ogt)F{V,t) 
as t 0. The Laplace transform version of Bohr’s formula then follows. 


(6.17) 

(6.18) 
□ 
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